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Abstract. We consider quantum-dynamical phenomena in the SU(2), S — 1/2 infinite-range quantum 
Ifeisenberg spin glass. For a fermionic generalization of the model we formulate generic dynamical self- 
consistency equations. Using the Popov-Fedotov trick to eliminate contributions of the non-magnetic 
fermionic states we study in particular the isotropic model variant on the spin space. Two complementary 
approximation schemes are applied: one restricts the quantum spin dynamics to a manageable number of 
Matsubara frequencies while the other employs an expansion in terms of the dynamical local spin suscep- 
tibility. We accurately determine the critical temperature Tc of the spin glass to paramagnet transition. 
We find that the dynamical correlations cause an increase of Tc by 2% compared to the result obtained in 
the spin-static approximation. The specific heat C{T) exhibits a pronounced cusp at Tc. Contradictory to 
other reports we do not observe a maximum in the C(T)-curve above Tc- 

PACS. 75.10.Nr Spin glass and other random models - 75.10.Jm Quantized spin models 



1 Introduction 

Theoretical models of spin glass systems, among which the 
best investigated is certainly the infinite-range Ising spin 
glass (SK model) J^, are well known for their conceptual 
difficulties due to peculiarities of the ordered phase. At- 
tempts of explicit solutions generally lead to severe numer- 
ical challenges. The infinite-range Heisenberg model, also 
known as the Quantum-SK model, is additionally com- 
plicated by the presence of quantum-dynamical correla- 
tions, even in the high-temperature phase. Realizations of 
Heisenberg spin glasses are discussed at length for exam- 
ple in the review by Binder and Young 0, emphasizing 
the role played by different types of anisotropy in contrast 
to the isotropic model. 

Theorists have been looking at the infinite-range Heisen- 
berg spin glass model from different angles: the quantum- 
dynamical self-consistency problem was first formulated 
by Bray and Moore in 1980 and the corresponding 
TAP-equations have been derived by Sommers |4l5j . Ef- 
fects of external fields and anisotropy were also investi- 
gated [H]. In these works explicit calculations relied on the 
spin-static approximation. Later, quantum-dynamical cor- 
relations were taken into account, for instance by means 
of Quantum Monte Carlo simulations of the paramagnetic 
phase [7] or by exact diagonalization of finite systems ^ . 
A different approach considers the SU(N)-generalization 
of the infinite range Heisenberg spin glass which can be 
solved exactly in the limit N ^ oo [911 1 1 . 

There are physically fascinating and further reaching 
questions concerning the interplay of frustrated magnetism 
and other interactions. The competition between the Hei- 
senberg spin glass interaction, Kondo-coupling, and trans- 



port, for instance, leads to non- Fermi liquid behavior jl2j . 
Despite the fact that such highly ambitious questions have 
been and continue to be addressed successfully, there re- 
main many not yet resolved basic problems even in the 
infinite-range Heisenberg spin glass model alone. 

This article is organized as follows. In sec.|21we intro- 
duce the fermionic model Hamiltonian and briefly sketch 
the dynamical spin glass decoupling procedure. We for- 
mulate general self-consistency equations for the fermionic 
SU(2), 5 = 1/2 infinite range spin glass model. Section 13 
concentrates on the isotropic Heisenberg spin glass on the 
spin space which is contained in the general model as a 
special case. After a short digression to the spin-static ap- 
proximation in sec. KHI we apply in sec. l3.2l the "dynamical 
approximation" which restricts the quantum spin dynam- 
ics to a feasible number of bosonic Matsubara frequencies. 
We accurately locate the paramagnet to spin glass phase 
transition in sec. 13. 2. II and present our results for the spe- 
cific heat in sec. 13.2.21 In order to support our results we 
additionally employ a perturbative expansion of the self- 
consistency equations in sec. 13.31 



2 Formulation of the general dynamical 
self-consistency problem 

2.1 Model definition and effective action 

We consider the grand-canonical Hamiltonian 
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where the latin indices label the N lattice sites and 
V = {x, y, z} denote the spatial directions. With the usual 
fermionic construction operators and a and the Pauli 
matrices cr^^ the spin 1/2 operators are represented by 



J2 



differs slightly from the one with couplings of the type 
Jij V ^^^^ been considered in other work 

3,4,5 6 7 8. In the paramagnetic phase, however, both 
model variants lead to identical effective actions and the 
results, particularly for the critical temperature, are thus 
directly comparable. 



where we dropped the conventional pre- factor H/2 for con- 
venience. We assume quenched disorder among the mag- 
netic couplings Jj^j — JJj according to the symmetric Gaus- 
sian distribution 



2ttJ, 



■ exp 



(2) 



The couplings in different spatial directions are completely 
uncorrelated. 

Although the Hamiltonian is defined on the Fock 
space it may readily be used to describe the corresponding 
model on the spin space, too. For this purpose one sim- 
ply has to choose the special imaginary and temperature- 
dependent chemical (Popov-Fedotov-) potential 



AiPF = -i7rT/2, 



(3) 



which effectuates exact cancellation of contributions of the 
superfluous non-magnetic (empty and doubly occupied) 
local states to the partition function |13| . 

We shall treat the general model (Jl within the frame- 
work of imaginary time Grassmann field theory, using the 
replica method along with standard decoupling techniques 
|3ll4j . Due to the infinite range of the magnetic interac- 
tions the system can be mapped easily onto a single-site 
problem by virtue of a site-global Hubbard-Stratonovich 
transformation. After a saddle point integration of the in- 
troduced decoupling fields the n-fold replicated and disor- 
der averaged partition function is given by a Grassmann 
path integral, 

[^"Idi.- JvWe^pi-A.s). (4) 
The resulting saddle point effective action reads 

i^f^iy ] 'Pt (5) 



where a and b are replica indices, ip' = if'^} is a Grass- 
mann spinor, S'^ = ^a^W denotes the Grassmann repre- 
sentation of a spin variable, and I2 is the 2x2 identity 
matrix. The r-integrations extend from to /? = 1 /T and 
the real functions Q'^l'^^i satisfy the saddle point condi- 
tions 



^ T — r' 



off 



(6) 



Note that the Hamiltonian ^ does not allow for qua- 
drupolar order, i.e. {S!^''^S'^f) tt = for v 7^ fi. Our model 



2.2 Self-consistency equations 

In the present article we consider properties of the model 
in the paramagnetic phase or in the spin glass phase close 
to the critical temperature such that a replica-symmetric 
treatment is sufficient. Hence, we introduce the compo- 
nents of a single time-independent spin glass order pa- 
rameter, 

Qtr' = 9- (7) 

The time dependence of the problem merely resides 
in the replica-diagonal spin-spin correlations Q^^'^/. We 
employ the Fourier decompositions of the time-dependent 
quantities. 



I — 00 

00 



(8) 
(9) 



where zi and tUm denote fermionic and bosonic Matsub- 
ara frequencies, respectively. The real Fourier coefhcients 
9m — = (^m) the Central quantities in our 

formulation of the theory. These parameters are closely 
related to the local dynamical spin susceptibility: 



(10) 



By means of a second decoupling step the interacting 
part of the effective action ij^l is rendered linear in the 
spin variables which permits the Gaussian integration of 
the Grassmann fields. For each spatial direction we in- 
troduce a replica-global decoupling field Zi, as well as a 
number of replica-local decoupling fields. The latter group 
into "static" fields ?/^,o and "dynamical" fields ?/^„j>i de- 
coupling spin-spin interactions which are diagonal and off- 
diagonal in the fermionic frequency space, respectively. 

This dynamical decoupling procedure leads to a system 
of non-interacting particles in the presence of an effective 
dynamical potential V. In the space of the fermionic fre- 
quencies zi , V takes the form of a non-diagonal Toeplitz- 
structured matrix defined by (for the sake of better read- 
ability we drop the replica indices from now) 



(11) 



v}^^ I' — I — m,m > 0, 



where the entries Vm themselves are 2x2 matrices in spin 
space. 
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Here the static effective magnetic fields are given by 

^^m=0 ^ K + Ju (V?^ + ^f^~q'v Vi^M^ , (13) 

whereas the dynamical effective fields are complex and 
read 

Kr>l = J^\l^{yt,m + iy^.m) ■ (14) 

As an auxiliary quantity we define the local one-particle 
Green's function in the presence of a particular configura- 
tion of these effective fields, 



V 



where 



(Gq ^)^,^ = {izi + fj.) (^r.zls 



(15) 



(16) 



is the inverse of the Green's function of the non-interacting 
system (i.e. Ji, =0). 

We condense the notation by the use of a space saving 
abbreviation for the Gaussian integral operator, 



G 



/(^)- (17) 



Furthermore, we introduce the following shorthand nota- 
tions for the multiple integrations over all occurring z-type 
and y-type decoupling fields: 



(18) 
(19) 




n n 



y--o m>i "'y-. 



Within the replica method, the disorder-averaged free 
energy per spin, /, is obtained by virtue of the relation 



which yields 



/3/ = -E 



lim 

n — >-oo 



- 1 



(20) 



suffices. 

By extremization of the free energy H21|) with respect 
to the parameters q^, and we finally derive the self- 
consistency equations 



(25) 



((Tr A;;r) (Tr A^„r) - Tr A'^.JA'^^r) (26) 
employing the auxiliary matrices 




iKi)i'i = cr^^Si,, 



(27) 



So far our formulation of the problem is rather general. 
With a suitable choice of the model parameters J^, h^, and 
^ eqs. (|ll|l - (|26|) can be used to investigate the Heisenberg- 
and XY spin glasses or the Ising spin glass in a transversal 
field, both on the Fock space as well as on the spin space. 
Note that the direction-dependent distribution ^ also 
allows for anisotropy. 



3 The isotropic Heisenberg spin glass on the 
spin space 

For the rest of this article we shall be concerned with 
the Heisenberg spin glass on the spin space. Therefore, 
the chemical potential will be fixed to ^ = /ipp given 
by eq. 0. We only consider the isotropic model in the 
sense that the distribution of the magnetic couplings Q 
are independent of the spatial direction, i.e. = J, and 
consequently so are the spin-spin correlations, i.e. q^, = q 
and q^ = q^. Furthermore, we restrict ourselves to the 
case of zero external fields. 

In anticipation of the following sections we state the 
important exact sum rule 



E 



m— — QO 



(28) 



where <P is given by 



which arises from eqs. Q and (|5J) at equal times due to 
the absence of non-magnetic local states in the model on 
the spin space. 



W. 



The weight function 

l¥ = det(r-i)/det(G74) 



(22) 3.1 The spin-static approximation revisited 



(23) 



results from the Gaussian integration of the Grassmann 
fields. For W to be finite and meaningful a regularization 
of the determinant in eq. H23|l is required. In this work the 
simple choice 

(G74)^,^=zza2 (24) 



As a starting point for the quantum-dynamical calcula- 
tions in sees. 13.21 and IBT^ and as a simple but instructive 
special case of the eqs. (|11I26|) we consider in this section 
the spin-static approximation 1^. 

The effective potential matrix (lll|l can be decomposed 
into a static and a dynamical part which are diagonal 
and non-diagonal in the fcrmionic frequency space, respec- 
tively: 

V = V,tat + Vdy„. (29) 
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According to the definition ^11^ Vstat is a block-diagonal 
matrix composed of the static fields whereas Vdyn 
comprises the complex dynamical fields (|14|l . 

In the spin-static approximation the time dependence 
of the saddle point ® is disregarded in the self-consistency 
problem. Consequently, the Fourier components q^^Q are 
neglected in the construction of the effective potential and 
hence Vdyn = 0. 

With the (frequency-diagonal) spin-static propagator 
matrix 

Tstat - (Gq-^ + Vstat)"' (30) 

the weight function can be calculated analytically 
resulting in 

W^stat = ^cosh(/3i/o), (31) 



where 



Ho 



(32) 



and the Hq are given by eq. Since Vdyn = the 

Gaussian integrations over the dynamical decoupling fields 
y^m>i '^^ ^'i^- H22|I " H26|I become trivial. If we restrict this 
discussion to the paramagnetic phase (where the spin glass 
order parameter vanishes, i.e. q — 0) eq. H22() evaluates to 



1 



' exp 1 y 



(33) 



with the abbreviation c = l3JyJ%. 

The dynamical approach made in this article facilitates 
the calculation of the dynamical saddle point components 
q^ even within the spin-static approximation. It can be 
seen easily that the first trace term in eq. H26|l vanishes 
for m 7^ for a frequency-diagonal matrix V . The second 
trace term, however, measures the overlap of two factors 
r that are displaced about m elements along the diagonal 
which yields 



3 ^ 



r A'^ r 

±1 AA^^^i stat/^— 7n' stat 



E 



2 {izi + n) {izi+jTL + m) 



2 tt2 
3^0 



7^2 



Performing the fermionic Matsubara sum and the angular 
integration of the j/i^^Q-fields (the Zj,-integrations become 
trivial in the paramagnetic phase) finally leads to 



3.2 Dynamical approximations 

In order to study the role played by quantum-dynamical 
correlations we adopt the method of dynamical approxi- 
mations that was introduced recently in the context of an 
itinerant spin glass model [TK] , 

In essence, this method systematically improves the 
spin-static approximation by successively taking into ac- 
count the dynamical contributions to the effective po- 
tential V Hll() . More precisely, in the so-called dynami- 
cal approximation of order M all Fourier components q^ 
with m = {0, .., M} are kept in the self-consistency struc- 
ture. The higher frequency components are set to zero 
in the construction of the effective potential. Thus, ujm 
plays the role of a cut-off frequency for the dynamical 
self-interaction. Technically, at order M the effective po- 
tential V is approximated by a band-diagonal matrix with 
band width 2(2M-h 1). 

The main benefit of this approximation scheme is that 
all Gaussian integrations over the dynamical fields J/j^^^jv/ 
become trivial which allows the numerical evaluation of 
the self-consistency equations H25|l and H26(l if the order 
M is small enough. 

The general strategy is to calculate a quantity within 
several dynamical approximations of increasing order M . 
Provided such a sequence of improved solutions reveals 
sufficient convergence properties the exact full dynamical 
result can be inferred by extrapolation to M — !■ cx). A sim- 
ilar approximation technique on the discretized imaginary 
time space has been used earlier in the context of the Ising 
spin glass in a transverse field |l()ll7j . 

In this article we present numerical solutions with M 
ranging from M ~ (spin-static approximation, see sec. 
13.1(1 to M = 4. Up to M = 2 the occurring integra- 
tions are performed by Gaussian quadrature. Due to the 
high dimensionality of the resulting integration problem 
for AI > 3 we applied a combination of Gaussian quadra- 
ture and a Monte Carlo method for some less important 
angular integrations. The latter is the origin of the statisti- 
cal errors in the numerical data presented in the following 
sections. 

The sum rule ((28|l was derived for the fully dynamical 
system. Nevertheless, it was verified to be fulfilled in the 
cases M — 1, 2 by high precision numerical calculations. 
Although we have yet been unable to prove it analytically 
for M > 0, we claim that eq. (|28|l holds exactly in any 
finite order of the dynamical approximation. 



1 



stat Jr 



cr sinh (cr) Smfi 



_|_ 7;-2j77,2 



cosh [cr] 



(34) 

For the static component this result constitutes a self- 
consistency equation with the explicit solution 



% 



-3 + + y/Q 30/32p -I- /34 j4 
6/32 J2 ■ 



(35) 



It is worth mentioning that the spin-static approximation 
in the shape of eq. ((34|l exactly fulfills the sum rule ((28|l . 



3.2.1 Solutions in the paramagnetic phase and 
determination of the critical temperature 

Our results for the zero-frequency part of the local sus- 
ceptibility xo H10|l obtained in the first three dynamical 
approximations are presented in fig. ^ We find that the 
quantum-dynamical corrections to xo relative to the spin- 
static approximation are quantitatively remarkably small 
(note the small vertical scale in fig. . This fact was al- 
ready pointed out in 7 . At high temperatures we observe 
quick convergence of this sequence of solutions meaning 



M. Bechmann, R. Oppermann: Dynamical solutions of a quantum Heisenberg spin glass model 



5 






0.6 0.8 1. 1.2 1.4 1.6 1.8 

TjJ 



Fig. 1. Paramagnetic solutions for the u)rn = part of the local 
susceptibility xo = Pqo in the dynamical approximations up to 
third order. Only deviations from the result in the spin-static 
approximation xo,o (eq. 13511 ') are shown. The dashed line rep- 
resents 1 — Jxo.o; according to eq. (I36II the intersection points 
mark the respective approximation to the critical temperature. 
To the right of the dashed line the shown paramagnetic solu- 
tions are unstable against spin glass order. 



0.4 - 


\ M = 




0.35- 


/\\ M = 2 




0.3- 






0.25- 






0.2- 







0.4 0.6 0.8 1. 1.2 1.4 1.6 1.8 



T/J 

Fig. 3. Specific heat according to eqs. 

lEHJ and lEZJ in the dy- 
namical approximations with M <2 (full lines). The cusp of 
the curves indicates the respective critical temperature Tc,m. 
Except for M = Q, in the spin glass phase the results are cor- 
rect close to Tc,M only (see text). For comparison the "conven- 
tional" spin-static approximation, where the frequency sum in 
eq. 13 711 is restricted to the m = term, is also shown (dashed 
line) . 




12 3 4 



M 

Fig. 2. Sequence of critical temperatures obtained within dy- 
namical approximations of orders M — {0, ..,4}. The open dia- 
mond and the dashed line represent the spin-static and the ex- 
trapolated full dynamical result, respectively. The arrows mark 
the critical temperatures obtained using the two lowest orders 
of the g„_^Q-expansionof the functional $ 12211 discussed in sec. 
13.31 Inset: Extrapolation to M = oo (open star) having regard 
to the 1/Af^-like convergence. The two thin lines indicate the 
statistical error. 



M — {0, ..,4} are presented in fig. From the structure 
of the self-consistency problem one expects a Af ^"^-like 
convergence of this sequence of Tc-approximants (see sec. 
lA.lfl . Extrapolation of the data to M — > oo, i.e. to the 
full quantum-dynamical solution, yields Tc — (0.58912 ± 
0.00015) J . This is an increase relative to the spin-static 
result Tcstat = 1/V3 J by about 2%. Our result can also 
be compared JS] to values for Tc obtained by means of 
quantum Monte Carlo simulations (Tc = 0.568 J) |Zj and 
exact diagonalization of finite systems (Tc « 0.52 J) 



3.2.2 The specific heat 

Using standard thermodynamic relations starting from 
the free energy H21|l we derive a useful expression for the 
internal energy per site, 

U = l^J' U - E qA , (37) 



that the quantum dynamics of the model is described vir- 
tually exactly by taking into account only the effects of 
very few Matsubara frequencies. Naturally, as the tem- 
perature decreases the number of the relevant frequencies 
increases. 

The equilibrium critical temperature of the paramag- 
net to spin glass phase transition can be determined by 
means of the simple relation 

Jxo{Tc) = l (36) 

which can be shown to hold in any order of the dynam- 
ical approximation by expansion of eq. (|25|l in powers 
of the order parameter q. Our solutions of eq. 1)36(1 with 



from which the specific heat 

CiT) = ^ (38) 

can be obtained by numerical evaluation of the temper- 
ature derivative. Within the dynamical approximation of 
order M the contributions to the frequency sum in eq. 
(|37|l with |m| > M were calculated non-self-consistently 
from eq. H26|l (or eq. (|34|l in the spin-static case M = 0). 

In order to investigate the behavior of the specific heat 
for T <Tc one needs solutions of eqs. H25|l and H26() in the 
spin glass phase which are easily found over the whole 
temperature range in the spin-static approximation. For 
M > 0, however, the full integration problem is hardly 
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feasible. Hence, we restrict ourselves to temperatures suf- 
ficiently close to Tc such that the self-consistency equa- 
tions can well be approximated by expansions in powers 
of Tc — T. From eq. (|25(l we obtain for the spin glass order 
parameter the linear expression 



with the slope 



1 
J 



dT 



(39) 



(40) 



Expansion of eq. H26|) in powers of q yields the simplified 
self-consistency equation 



— Rr, 



T < T„ 



(41) 



where R„i symbolizes the right hand side of eq. (|26|l and 
the Cm are well defined expansion coefficients that can be 
calculated numerically at T — Tc. By using the solutions 
of eqs. and lHT|l in eq. l|57|l we obtain curves for the 
specific heat in the ordered phase that are correct at linear 
order of Tc — T. Note in particular that effects of Parisi 
replica symmetry breaking will change the results only in 
higher orders of T = Tc. 

The resulting specific heat approximants with M = 
{0, ..,2} are shown in fig. |3| Due to the apparent quick 
convergence of this sequence of solutions we may safely 
draw qualitative conclusions for the limit Af — > oo. In 
the paramagnetic phase C(T) monotonically increases as 
the temperature is lowered. Contrary to what was previ- 
ously reported by other authors 8 in our results there 
is no indication of a broad maximum in the full dynami- 
cal C(T)-curve above Tc. Merely the "conventional" spin- 
static approximation which neglects the quantum dynam- 
ics altogether and omits all m ^ terms in the internal 
energy formula l|37|) has such a feature (fig.El)- 



3.3 Perturbative expansion of the self-consistency 
equations in powers of the dynamical parameters q^^q 

To consolidate our findings of sec. 13. 21 we additionally ap- 
plied a completely different approximation technique to 
the self-consistency problem. Viewing the spin-static the- 
ory of sec. 13. li as a starting point we perturbatively expand 
the functional (I22|l in powers of the dynamical saddle 
point components q^^o- Such an expansion is certainly 
justified at high temperatures where q^-^Q ^ 9o ^ 1- 
With eq. (|21|l the (7,„_^Q-expansionyields an approximate 
free energy which has to be extremized with respect to the 
parameters q^. Thus, we derive simplified self-consistency 
equations fsee lA.ip that contain only integrations over the 
static decoupling fields yi,^o and are therefore easily solved 
numerically. 

The resulting solutions for the zero frequency suscepti- 
bility and the specific heat at the two lowest orders of the 
g^-^Q-expansionO (q^) and O {q^Qm') shown in fig.^ 
and fig. 13 respectively. At high temperatures the previous 



I 0.008 

o" 

X 

•"5 0.006 
I 

§ 0.004 



C(9m3m') 



T/J 

Fig. 4. Paramagnetic solutions for xo ~ P<io (see fig. in 
the first two orders of the g^^o'^^PS-nsionand in the third or- 
der dynamical approximation (plot symbols) for comparison. 
Intersections with the dashed line mark the respective approx- 
imate critical temperatures which are indicated in fig. |5| 



C „. 




Fig. 5. Specific heat in the first two orders of the g„-expansion 
and in the second order dynamical approximation (dashed line) 
for comparison (see fig. O . Only solutions in the paramagnetic 
phase are shown. 



.J 




Fig. 6. Check of the sum rule 12811 for the two lowest orders of 
the g^-^o'Sxps-nsionindicating the quality of these approxima- 
tions. Note that the higher order O {qmQm.') clearly does better 
for all temperatures. 
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results of the dynamical approximation discussed in sec. 
13. 21 are reproduced very accurately. However, for T/ J < 1 
the results obtained within the two different approxima- 
tion schemes clearly differ from each other and apparently 
higher orders of the g^^_^g-expansionbecome important. 

From the numerical data presented in this section it is 
not obvious that the two sequences of approximations, the 
dynamical approximation with increasing M on one hand, 
and the increasing orders of the g^__iQ-expansionon the 
other hand, will finally converge to the same full dynam- 
ical solution as required. To prevent miss-interpretations 
we check in fig. how accurately the sum rule is 
obeyed by the g^_^Q-expansionat the two lowest orders. 
At high temperatures the sum rule is fulfilled almost ex- 
actly at order O [q^q^y^i) reflecting the good quahty of the 
approximation. However, we observe a substantial viola- 
tion of the sum rule for T / J < 1 providing clear evidence 
that in this temperature regime this sequence of solutions 
is not well converged yet at the highest order considered 
in this article. 

Due to the small number of available orders we can not 
reliably extrapolate the results of the 5^_^Q-expansionto 
the full dynamical quantities. It can be seen in figs. El 
and ^ however, that our results for the zero- frequency 
susceptibility presented in sec. 13.2.11 and particularly our 
statement for Tc are strongly supported by the present 
perturbative calculation. The specific heat, on the other 
hand, directly depends on the dynamical parameters <7„-^o 
and is therefore more sensible to the failure of the approx- 
imation to fulfill the sum rule (|28|l . We can nevertheless 
draw qualitative conclusions. The maximum obtained in 
the "conventional" spin-static approximation (see fig. ^ 
being the zeroth order of the (7,„_^Q-expansionis already 
very weak at order O {qm)i it is not present any more 
at order ©(g^g^,). Hence, consistently with sec. 13.2.21 
we claim that there is no maximum in the full dynamical 
C(T)-curve above T^. 



4 Summary and conclusion 

In this article we applied standard techniques for many 
body systems to the fermionic SU(2), S — 1/2 infinite 
range spin glass model and formulated the general dy- 
namical self-consistency problem for the spin-spin corre- 
lations in secEl Using the Popov- Fedotov potential Q to 
eliminate contributions of the non-magnetic local states 
we studied the isotropic Heisenberg spin glass on the spin 
space in sec. 13 The results for the corresponding fermionic 
model on the Fock space will be published elsewhere. 

In order to solve the highly coupled self-consistency 
equations and particularly to make the high-dimensional 
integration problem in eqs. H25|l and (|26l) feasible we used 
two different systematic approximation schemes. The dy- 
namical approximation of order M, on one hand, describes 
the quantum dynamics with a limited number of Matsub- 
ara frequencies but the corresponding saddle point com- 
ponents 'Z|„j|<jv/ are dealt with exactly fsec. l3.2|) . The per- 
turbative (7j^_^Q-expansion, on the other hand, takes into 




Fig. 7. Diagrams contributing to the first two orders of the 
perturbative g^-^g'^^PS-nsionof the functional <1> defined by eqs. 
12211 and 1421 1. Straight lines represent the spin-static propaga- 
tor wavy lines symbolize contractions of the dynamical 
efi'ective fields dJ. 



account all frequencies but only a few powers of g„j-^Q 
(sec. 13.311 . In this sense the two approximation schemes 
are complementary to each other. 

Both approaches yield a consistent picture for the zero 
frequency local spin susceptibility in the paramagnetic 
phase (figs, n and 01). By extrapolation of the results in 
the dynamical approximation of the orders M — {0, ..,4} 
to M ^ oo (fig. O we obtained the full dynamical criti- 
cal temperature = (0.58912±0.00015) J which is about 
2% higher than the value in the spin-static approximation. 
We also presented results for the specific heat C{T) (figs. 
and [SJ . In the framework of the dynamical approxima- 
tion we perturbatively extended the calculations to the 
spin glass phase. We found a pronounced non-analyticity 
of C{T) at Tc- We can not confirm the observation of a 
broad maximum in the C(T)-curve above Tc which was 
reported by other authors (HI- In this point our results 
rather resemble those for a SU(iV)-generalization of the 
quantum Heisenberg spin glass in the limit N oo and 
for larger S [TTITTI . 

The methods discussed in this article are particularly 
useful to qualitatively and quantitatively describe the high 
temperature phases of disordered quantum systems. There 
are many open issues that can thus be addressed, e.g. real- 
frequency response functions, the behavior in a magnetic 
field, or questions concerning anisotropy. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft under research project Op28/5-2 and by the 
SFB410. One of us (M. B.) also wishes to acknowledge the 
scholarship granted by the University of Wiirzburg. 



Appendix 

A.l Details of the perturbative (2',jj_^Q-expansion 

As a consequence of the decomposition (|29ll of the effective 
potential (|ll|l the weight function (|23() can be written as 

W = Wstat exp Tr In (1 + T.tat Vdyn) , (42) 

where the spin-static quantities Fstat and Wstat are defined 
by eqs. iPHjl and respectively. 
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In the high temperature limit the dynamical saddle 
point components vanish Uke (j„-^g ~ whereas the 

zero frequency component reaches unity, % ^ 1. Hence, 
at high temperatures the logarithm in eq. 142|l can be 
expanded in powers of the matrix Vdyn- The trace in eq. 
(|42(l and the Gaussian integrations over the dynamical 
decouphng fields y^„>i in eq. (|^ yield a representation 
of the functional 



<?1 + <?2 + 0{q^q^,q„ 



(43) 



in terms of contracted diagrams. 

Figure displays all occurring diagrams up to the sec- 
ond non-trivial order. Evaluation of these diagrams yields 
the final results 



^ cr sinh (cr) 



m—1 



(44) 



and 



OO 

^rnQm' (w Cr, (tTTTi) 



'>T1 



cr, (7rm)^ , {Trm'Y 



(45) 



with c = (3J yjq^ and the functions 

secha; / Ssx^ — s^.t^ + 2s'^ 



(a;2 + y (2;2 s) (a;^ + Asx) 



sinh (2a;) 



4 , , lla:^ 
— cosh (2a;) H 



3s -4 



(46) 
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sinh a; ( 

' ' {x^ + s) {x^ + t) \ cosh X sinh a; x"^ -\- s 

2tx X (5a;2 + s + t){3 + t) 



a;2 -I- t (a;2 + .s + t) - Ast 



(47) 



Expressions (|43|I - H47|I correctly reproduce the high tem- 
perature expansions of all quantities including order O (/3^) . 

In the context of the dynamical approximations dis- 
cussed in sec. 13.21 particularly for the extrapolation to the 
full dynamical result of a quantity it is important to know 
how this quantity varies with the order M for M oo. It 
is sufficient to restrict the discussion to (|44ll . We write 

= <?i,j\/ + <^i,M, (48) 

where for ^i^m the sum in eq. (|44|l is restricted to m = 
{1, .., M} and 'Pi^m contains the remaining high- frequency 
terms with m > M. Since ^ m^^ foj- large m (see eq. 
(EH)) it is evident that Si^m ~ M'^ for large M. Thus, 
the functional and consequently all derived quantities, 
e.g. the critical temperature, converge like M~^. 



